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1.0  INTRODUCTION 


The  TRAM  program  provides  a  Monte  Carlo  simulation  to  determine  a  cumula¬ 
tive  hit  probability  for  an  ASW  weapon  delivery  system  with  an  evasive  sub¬ 
marine  as  target.  The  target  evasive  maneuver  is  a  course  change  with  a 
simultaneous  speed  change  order  ( new  thrust  level ) .  By  user  choice ,  the 
maneuver  parameters  may  be  fixed  or  may  be  random  within  prescribed  statisti¬ 
cal  limits.  The  program  also  takes  into  account  random  errors  of  the  weapon 
delivery  system.  These  include  errors  in  target  localization,  target  motion 
analysis  (TMA) ,  and  the  delivery  error  of  the  system. 

Hit  probability  for  a  given  weapon  delivery  point  is  determined  by  a 
user-provided  rectangular  grid  of  hit  probabilities.  The  grid  structure 
provides  rectangular  areas  distributed  about  the  target  with  a  prescribed  hit 
probability  associated  with  weapon  delivery  in  each  area. 

The  program  is  written  in  Sperry  Uni vac  1100  Series  FORTRAN  (ASCII)  — 
Level  9R1 .  See  reference  1 .  A  complete  source  program  listing  is  contained 
in  Appendix  A. 

Logical  simplicity  of  the  source  code  was  emphasized  to  make  the  program 
and  its  logical  flow  more  easily  understood  by  the  user.  In  some  instances, 
this  resulted  in  deliberately  redundant  or  inefficient  source  code. 


2.0  MATHEMATICAL  DESCRIPTION 

Figure  1  presents  a  representative  illustration  of  target  motion  and 
weapon  delivery  system  geometry.  The  geometry  is  two-dimensional »  target 
depth  is  assumed  constant. 

At  time  t  *  0  the  target  is  located  at  position  T(0),  located  at  the 

origin  of  the  xy  coordinate  frame,  with  an  initial  velocity  v  along  the 

o 

positive  x-axis.  At  some  subsequent  time  t  ■  T«,  the  target  initiates  a 

(1]  Sperry  Rand  Corporation.  Sperry  Univac  Series  1100  FORTRAN  (ASCII)  - 
Level  9R1,  Programmer  Reference.  1979. 
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Figure  1 .  Illustrative  run  geometry. 


maneuver  which  consists  of  a  course  change,  or  an  instantaneous  change  in 
thrust  level  (change  in  ordered  speed),  or  both,  or  neither. 

The  time  t  =  T1  may  be  prescribed  by  the  user  or,  otherwise,  will  be 
chosen  by  the  program  from  a  random  distribution  of  times  uniformly  distri¬ 
buted  over  the  interval  from  t*0tot  =  Tg  +  e^,  where  Tg  is  the  estimated 
blind  time  of  the  weapon  delivery  system,  ie,  the  time  between  last  observed 
target  position  and  weapon  delivery,  and  e,^  is  the  error  in  the  estimate  of 

T  .  e  is  a  randan  variable  chosen  by  the  program  from  a  normal  distribu- 
B  To 

tion,  with  zero  mean  and  a  user-specified  standard  deviation. 
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The  thrust  level  change  at  t  ■  T^  Is  specified  by  the  user.  The  course 
change  magnitude  may  also  be  specified  by  the  user  or,  otherwise,  is  chosen  by 
the  program  from  a  uniform  distribution  between  limits  provided  by  the  user. 
The  direction  of  the  course  change  is  randomly  chosen  by  the  program,  with 
either  direction  equally  probable. 

The  maneuver  event  itself  may  be  treated  as  a  random  variable  by  the 
program.  That  is,  while  the  user  may  prescribe  that  a  maneuver  does  occur,  he 
may,  alternatively,  provide  only  a  probability  that  a  maneuver  occurs.  In  the 
latter  case,  the  program  will  determine  the  occurrence  of  the  maneuver  by 
comparing  a  random  number,  drawn  from  a  uniform  distribution  over  the  interval 
zero  to  one,  with  the  prescribed  probability. 

The  time  t  ■  is  the  time  at  which  the  course  change,  if  any,  is  com¬ 
pleted.  Thereafter,  the  target  will  proceed  on  a  steady  course  with  constant 

thrust  level  until  the  end  of  the  run  at  time  t  *  T  +  e _ . 

B  TB 

Target  dynamics  during  the  maneuver  are  determined  from  the  two-degree- 
of-freedom  equations  described  in  Appendix  B. 

At  time  t  -  0,  the  weapon  delivery  system  assumes  the  target  to  be  at 
position  P(0).  e^  and  e^  are  the  x  and  y  target  localization  errors.  These 
errors  are  chosen  by  the  program  from  a  zero  mean  normal  distribution  with  a 
user-specified  standard  deviation. 


The  weapon  delivery  system  assumes  that  the  target  maintains  a  steady 

course  e^  and  steady  speed  vq  +  e^  from  time  t  ■  0  until  t  m  Tg.  The  target 

course  estimate  error  e  is  chosen  by  the  program  from  a  zero  mean  normal 

9 

distribution  with  a  standard  deviation  specified  by  the  user.  The  target 
speed  estimate  error  e^  is  chosen  from  another  zero  mean  normal  distribution 
with  a  user  specified  standard  deviation. 

At  time  t  -  T  ,  the  weapon  delivery  system's  estimated  (predicted)  target 
B 

position  is  at  position  P(T  ).  The  aim point  AP  is  specified  by  the  user  in 
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terns  of  range  ( R^)  and  bearing  (♦ft)  offsets  relative  to  the  predicted  target 
position  at  t  =  Tg. 

The  weapon  delivery  point  is  at  DP.  The  delivery  point  is  offset  from 
the  aimpoint  AP  by  the  weapon  delivery  errors  e^  and  e^.  Both  are  chosen  by 
the  program  from  a  zero  mean  normal  distribution  with  a  user-specified  stan¬ 
dard  deviation. 

The  probability  of  hit  for  a  given  weapon  delivery  is  determined  by  the 
location  of  the  weapon  delivery  point  in  a  user-defined  hit  probability  grid. 
The  rectangular  grid  is  formed  by  a  series  of  straight  lines  parallel  to  the 
target's  course  at  weapon  delivery  and  by  another  series  of  straight  lines 
perpendicular  to  the  target  course.  An  illustrative  example  is  shown  on 
Figure  2. 

The  target  is  located  at  the  origin  (x  =  0,  y  =  0)  and  the  target  course 
is  along  the  positive  x-axis.  The  grid  lines  are  specified  by  the  user  and 
may  be  located  anywhere  relative  to  the  target.  There  is  no  requirement  for 
symmetry  about  any  axis. 

The  user  provides  a  probability  of  hit  for  each  rectangle  formed  by  the 
grid.  The  left-hand  and  lower  boundaries  of  each  rectangle  are  assumed  to 
belong  to  that  rectangle.  Thus,  for  example,  if  the  weapon  delivery  point  was 
located  anywhere  in  the  rectangle  defined  by  500  ^  x  <  1000  yards  and  -900  ^  y 
<  -500  yards,  the  hit  probability  for  that  weapon  delivery  would  be  0.8. 

An  overall,  cumulative  hit  probability  is  computed  from  a  set  of  hit 
probabilities  from  individual  "runs",  where  a  run  is  defined  as  one  target 
trajectory  and  associated  weapon  deliveries.  A  "set"  of  runs  is  defined  as  a 
collection  of  runs  made  under  conditions  that  differ,  one  run  from  another, 
only  in  a  stochastic  sense.  This  structure  comprises  a  Monte  Carlo 
simulation. 
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The  user  specifies  the  number  of  runs  in  a  set  and  may  also  specify  more 
than  one  set  of  runs  for  a  given  program  execution.  He  may  also  specify  any 
number  of  weapon  deliveries  for  each  run  (salvo  launch). 

For  a  run  with  salvo  launch  of  N  weapons,  the  hit  probability  for  the  run 
is  computed  according  to  the  equation: 


P  =  P  +P.-PP, 
n  n  n-1  n  n-1 


n  *  1,  2| . . .H 


where  P  is  the  probability  that  at  least  one  of  the  first  n  weapons  will  hit 

'  n  th 

and  P^  is  the  hit  probability  for  the  n  weapon. 

The  cumulative  hit  probability  (Monte  Carlo  probability)  for  a  set  of  N 
runs  is  computed  according  to  the  equation: 


i  E 


where  P^  is  the  hit  probability  for  the  k  run  in  the  set. 

The  random  variables  for  each  run  are  provided  by  two  random  number  gen¬ 
erator  subroutines:  UNFRM  and  GAUSS ,  Both  generators  run  continuously  (with¬ 
out  reset)  through  a  set  of  runs,  so  generally  different  sequences  of  numbers 
are  provided  for  different  runs  in  a  set.  A  seed  number  for  the  generators  is 
provided  by  the  user  at  the  beginning  of  each  set  of  runs. 

UNFRM  generates  a  sequence  of  real  numbers  uniformly  distributed  over  the 
interval  (0,1].  The  algorithm  used  is  a  multiplicative  type  of  the  form: 


x  «  (185363)  x 
n+1  n 


(mod  235) 


The  interval  of  the  distribution  is  changed  to  (0,1]  by  discarding  all 


27  -27 

x^  >  2  and  multiplying  the  others  by  2 


The  2  factor  is  a  consequence 


of  the  integer-to-real  conversion  process  used  with  the  Sperry  Univac  1100 
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Series  FORTRAN  (ASCII)  compiler.  See  references  1  and  2  for  additional 
details . 

The  subroutine  GAUSS  generates  a  sequence  of  normally  distributed  random 
numbers.  The  algorithm  is  based  upon  the  Central  Limit  theorem  and  is  of  the 
form; 


where  x  is  the  normally  distributed  random  variable,  the  y^  are  random  numbers 

uniformly  distributed  over  the  interval  [0,1]  (provided  by  subroutine  UNFRM) , 

o  is  the  standard  deviation  of  the  distribution,  and  x  is  the  mean,  x  and  o 
x  x 

are  specified  in  the  program  call  for  the  subroutine.  See  reference  2  for 
additional  details. 

Source  code  listings  of  UNFRM  and  GAUSS  are  included  in  Appendix  A. 

3.0  INPUT  DATA  REQUIREMENTS 

For  the  convenience  of  the  user  in  formatting  input  data,  all  input  data 
variable  names  are  defined  with  a  NAMELIST  statement.  See  reference  1  for 
details  of  NAMELIST  input. 

The  program  defines  three  NAMELIST  names:  NL1,  NL2 ,  and  NL3.  Every  input 
data  variable  is  assigned  to  one  of  the  three.  Upon  initiation  of  program 
execution,  the  program  searches  for  and  reads-in  NL1  data.  Upon  initiation  of 
each  set  of  runs,  the  program  searches  for  and  reads-in  NL2  and  NL3  data.  For 
the  first  set  of  runs,  all  NL2  and  NL3  data  variables  must  be  explicitly  de¬ 
fined  on  an  input  record  (card).  For  subsequent  sets,  only  those  input  vari¬ 
ables  whose  values  change  from  the  preceding  set  need  to  be  redefined.  If 
none  of  the  variables  in  NL2  or  NL3  change  value  from  the  preceding  set,  an 
input  card  with  a  blank  data  field  must  nevertheless  be  provided  for  that 


[2]  Hanning,  R.  W.  Numerical  Methods  for  Scientists  and  Engineers.  McGraw- 
Hill,  1973. 
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NAMELIST  name,  since  the  program  will  search  for  both  NL2  and  NL3  input  cards 
before  executing  a  set  of  runs. 

Appendix  C  contains  a  listing  of  all  input  data  variable  names  and  a 
description  of  the  datum  each  represents.  Appendix  D  contains  an  example  of 
input  data  cards  to  illustrate  input  data  requirements  and  format. 

4.0  SUBROUTINE  REQUIREMENTS 

The  main  program  calls  two  subroutines:  UNFRM  and  GAUSS.  Subroutine 
UNFRM  (NO,  Nl,  X)  returns  to  the  calling  program  a  real  number  X,  which  rep¬ 
resents  one  pseudorandom  sample  from  a  uniform  distribution  over  the  interval 

[0,1],  and  an  integer  number  Nl,  which  represents  one  pseudo-random  sample 

27 

from  a  set  of  integers  uniformly  distributed  over  the  interval  [0,2  ].  Sub¬ 

routine  GAUSS  (NO,  XMEAN ,  XSIGMA,  Nl ,  X)  returns  to  the  calling  program  a  real 
number  X,  which  represents  one  pseudorandom  sample  from  a  normal  distribution 
with  mean  XMEAN  and  standard  deviation  XSIGMA,  and  an  integer  number  Nl ,  which 

represents  a  pseudorandom  sample  from  a  set  of  integers  uniformly  distributed 

27 

over  the  interval  [0,2  ] . 

For  both  UNFRM  and  GAUSS,  the  integer  NO  is  provided  by  the  calling 
program.  For  the  first  call  for  either  UNFRM  or  GAUSS,  NO  may  be  any  odd 
integer.  If  either  UNFRM  or  GAUSS  has  been  called  prior  to  a  current  call  for 
either  subroutine,  then  NO  should  be  the  value  of  Nl  returned  from  that  call, 
for  either  subroutine,  which  immediately  precedes  the  current  call. 


'  i 
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Source  code  listings  of  UNF»!  and  GAUSS  are  included  in  Appendix  A.  A 
mathematical  description  was  provided  in  Section  2.0.  Note  that  subroutine 
GAUSS  calls  the  subroutine  UNFRM. 


5.0  OUTPUT  DATA 

Appendix  E  contains  the  printout  for  program  execution  with  the  input 
data  illustrated  in  Appendix  D. 
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6.0  TOP  LEVEL  FLOWCHART 


The  top-level  program  flowchart  is  shown  in  Figure  3.  At  the  level 
shown,  program  execution  is  essentially  controlled  by  three  nested  do-loops. 

The  outermost  loop  (1=1,  NSET)  cycles  once  for  each  set  of  runs,  while 
the  second  loop  (J=1,  NRUN)  cycles  once  for  each  ran  within  a  set.  The 
innermost  loop  (K=1,  NSHOT)  cycles  once  for  each  weapon  shot  of  a  run.  With 
one  exception,  lower  level  flowcharts  are  not  necessary,  since  program  logic 
is  obvious.  A  second-level  flowchart  for  the  one  exception.  Target  Course  and 
Position  at  Weapon  Delivery,  is  presented  in  the  next  section. 

7.0  TARGET  COURSE  AND  POSITION  AT  WEAPON  DELIVERY  -  LEVEL  II  FLOWCHART 

Figure  4  is  a  level  II  flowchart  for  target  course  and  position  at  weapon 
delivery.  The  logical  flow  of  program  execution  can  be  discerned  from  the 
flowchart  and  the  comments  included  in  the  source  code  listing  (Appendix  A). 

8.0  INTERNAL  VARIABLE  AND  FUNCTION  NAMES 

Appendix  F  contains  a  listing  of  all  internal  variable  and  function  names 
used  in  the  program  along  with  a  d^inition  of  each. 
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Figure  4.  Target  course  and  position  at  weapon  delivery  -  Level  11  Flow  chart 
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standaRO  deviation  of  circular  normal  Error 

38. 

c 

of  WEaRON  DELIVERY 

39. 

c 

NX 

•  •  m  • 

number  OF  MIT  PROBABILITY  GRID  LINES 

90. 

c 

PERPENDICULAR  TO  TARGET  CENTERLINE  125  mAX) 

91  . 

c 

NY 

•  ••  • 

NUMBER  OF  HIT  PROBABILITY  GRID  LINES  PArALLEL 

92. 

c 

TO  TARGET  CENTERLINE  (25  maxi 

93. 

c 

NRAN 

•  m  •  9R 

ARBITRARY  INTEGER  SEED  F OR  RANDOM  NUmBEr 

9  9  * 

c 

GENERATOR 

95. 

c 

GP|DX( 

1  )  YOS 

ARRAY  (25  MAXI  OF  X  COORDINATES  Of  MiT 

96. 

c 

PROBABILITY  GRID  LINES.  POSITIVE  AHEAD, 

97. 

c 

NEGATIVE  astern.  ARRAY  IS  STRUCTURED  IN 

99. 

c 

ASCENDING  ORDER  *  GRIDXII)  »LT.  GR|DX(l*ll, 

9  9. 

c 

GR  !  t>  Y  (  1  )  YOS 

ARRAY  (25  MAX)  CF  V  COORDINATES  OF  H|T 

so. 

c 

probability  grio  lines.  POSITIVE  TO  PORt, 

si  . 

c 

negative  to  stbo,  array  is  structured  In 

52. 

c 

ASCENDING  ORDER  -  GR  I  OY  (  |  I  ,  LT  •  GRIDYlUt). 

53. 

c 

PH(  1 

ARRAY  (29  BY  29  MAX)  OF  HIT  PROBABILITIES. 

59. 

c 

PH(I.J)  IS  THE  PROBABILITY  ASSOCIATED  WjTH  THE 

55. 

c 

hit  probability  grid  rectangle  defined  rt 

-2 


c 

c 

c 


IlGRIOXIIl  ,LE.  X  I  ,ANt.  (  G  H  I  0  X  I  |  ♦  J  )  ,Gt.  x  I  I 

,  AND. 

IlGRlDYtJl  .LE  •  Y>  ,ANt.  tGRlOYlJM)  ,GT.  Y>> 


54 

57 

s* 

59 

60 
6) 

42 

43 
49 
45 

44 
47 
4ft 

69 

70 
7  J 

72 

73 
79 
75 

74 

77 

78 

79 

80 
81 
82 
83 
89 

85 

86 

87 

88 

89 

90 
9  I 
9  2 
93 
99 

95 

96 

97 

98 

99 
100 
I  C  1 
102 
103 
109 

105 

106 
107 
ICR 

109 

110 
I  1  1 
I  12 


C  For  THf  FIRST  SFT,  ALL  INPUT  DATA  MUST  PE  PROVIDED.  FOR  SUBSEQUENT 

C  SfTS,  OILY  data  motificai ions  are  required,  nset  is  read  only  ONCE 

c  at  twe  start  of  the  program,  all  other  data  are  read  at  start  of 

C  EACH  set,  namelist  INPUT  is  used. 

f  SPECIFICATION  STATEMENTS 

DIMENSION  GR | OX  I  25  I  ,  GR 1 0  Y ( 25  I  |Pn(29,?9I 

namelist  /nh/nset 

X  /NL2/NRUN  jNSMOT .PTM.TSl  ,TS2, CMAX.CM1N.tr, TCL.TR2.Tb. APR, 

X  APA.SIGTB.SIGTC.SIGTS.SIGL ,SJGD,NX,NV,NRAN 

x  /NL3/GRlDX,GRir)Y,PH 

0  22  FoRMTflHj  ,  99HCUMUL  AT !  VE  HIT  PROBABILITY  FOR  ALL  RUNS  THIS  SETS. 

X  2X.6HPCH  ■  ,F7>//t*,30HlNPliT  DATA  FOR  THE  SET  FOlLOwII 

026  format i i h  ,i9hgrioxii)  i 1-1 , , «2,2h> :/ioifs. 1 :2*  1 1 

027  F  C  R  M  A  T (  I H  ,19HGRI0Y(J|  ( J ■ 1 .  .12, ?H)  5 / 1 0 ( F  8  .  |  ; 2  X  )  | 

028  FCRmATIIH  ,13HPH(I,J)  C I « | , , 1 2 , 7N I , I J* I , , 1 2 , 2H 1 ; / | 0  I F8 • 9 » 2*  I  > 

c  define  inverse  hyperbolic  functions 

ASU)H(  X  1  *ALOG<  X ♦SORT  I  **X*  |  .0)  1 
ACOSH| X ) ■ ALOG ( X*SQRT ( X*X- 1 .0  >  > 

AT  ANN  I  X  I »0.5« ALOGI  (  I .o*X 1 / I I.0-X1  I 
AC0tH(XI«O.5»AL0G(  ( X  ♦  i ,p)/(X»l .01  1 

c  initiate  loop  for  sets  of  runs 
readis.nli  i 

DC  019  !»l  , NSE T 

c  read  input  data  for  set 

Re  ad  1 5 .ml 2 1 

READI5.nl 31 

C  INITIALIZE  DATA  FOP  this  SFT 

Dun  1*0.562962963 
TSI I *DUM I • T  S 1 
TS2 I “OUM | • TS2 
SlGrSl»PUMl«S|GTS 
TCL I »T Cl/3.0 
OUM, «o.OI 795329 
APAJ»DUMi*APA 
SlGTCI«OUM|«SIGTC 
NpAn 1 «?*NR AN* | 

PCH.o.O 

FnRuN«FLO A  T ( NRUN ) 

c  initiate  loop  for  runs  in  this  set 

DO  017  J  « I »  N  R  U  N 

C  TIME  TO  WEAPON  OELIvERY  for  this  RUN 

Call  GA<  SSINRAN1 ,0.0.5IGTB,NRAN|  ,T| ) 

T I *TP*T I 

C  COMpUTFD  TARGET  POSITION  AT  WEAPON  DELIVERY 

CALL  GAUSS  I  NR  AN  I  , 0 . 0 . S I GL , NR  AN |  ,  X TC I 
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1 3 
19 

15 
t  6 
1  7 

ie 

1  9 

70 

71 

72 
23 
29 

25 

26 
27 
2fl 

29 

30 

31 

32 

33 
39 

35 

36 

37 
30 
39 
*•0 
9  I 
97 
93 
99 

95 

96 

97 
90 
99 

50 

51 

52 

53 
59 

55 

56 

57 
50 

59 

60 
61 
62 
63 
69 

65 

66 
67 
60 
69 


call  fiA"ssiNP*Ni .o.c.sigl.nram, ,»tci 
call  gauss  i  nr  «m .o.o.sigtsI .npani ,dumi  i 

OUHj«»  T *•  l  I  ♦  0  U  M  |  I  •  T 0 

CALL  GAl'SSlNPAM  .o.o.SlGTCl  ,NR*M! 

XTC»XTf*DU*i»C0S(0UM21 
VTC»YTe-DLN j  «S JN1DUM2 | 

C  AJMpOINT 

0UM2«0tlM2*APA  t 
XAP»XTC*APR»C05I0UM2) 

YAP»VTC-APR*S1M(0U»'2  1 

C  TARGET  COl'PSE  ANO  POSITION  AT  WFAPON  DELIVERY 

C  IS  TMfRE  A  TARGET  MANEUVER 

call  unfkm<nrani .nrani ,rn> 

in | I PTN.LT.O.O) .OR. IPN.LT.PTH) 1 . AND. « ( CM A  X . NE • 0 . 0 1 #0R. 

X  <  TS2.»  E.TSj  )  |  1  GO  TO  CO | 

C  TARGET  COURSE  AMO  POSITION  AT  WEAPON  DELIVERY.  NO  MANEUVE*. 

TCWO-O.P 
X  T*  TS  1  T • T I 

Y  T  »0 • 0 

GO  TO  002 

C  TARGET  maneuver,  time  maneuver  starts. 

roi  IF  1PTM.LT. 0.0)  GO  TO  n I  2 
call  unerminrani ,nrani ,rni 

T  M  ■  R  M  *  T  I 
GO  TO  013 

C I  2  Tm«-| ,o«PTM«T| 

013  CONTINUE 

c  target  position  at  start  oE  maneuver 

XT*TS| i«tm 

Y  T  *0  •  0 

C  OCES  MANEUVER  INCLUDE  COURSE  CHANGE 

Call  unfrmjnrani .nrani ,rni 

CCM«CMIN*RN*(CMAX-CM1*II 
IFICCn.NE.O.Cl  GO  TO  003 

C  TARGET  COURSE  AMO  POSITION  AT  WEAPON  DELIVERY.  SPFED  CHANGE  ONLY, 

TCWD-O.u 
FNUO-TSI /TS2 
TAU.TS2 I • 1 T l-TM 1 /TCLI 
IF( FNdo.GT.  |  .o  |  GO  TO  009 

C  TARGET  ACCELERATES 

PUMi»ATANH(FNUOI 

XT»XT»TCLI»ALOGICOSH(TAU*DUMi I/COSHIDUM, ) ) 

GO  TO  002 

c  target  decelerates 

009  OuM j BACDTHl FNUO 1 

XT»XT*TCL I • A  LOG ( S(NH( TAU*OUmI  l/SINHCOUHl  I  1 
GO  TO  01)2 
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TARGET  COURSE  AND  POSITION  AT  WEAPON  DELIVERY.  COURSE  CM*NGr 
0"3  Ot)M|-|  .0*Tf?2/<  TR«TR| 

CR»TCLI/OUN| 

MAGNITUDE  AND  DIRECTION  OF  COURSE  CHANGE 
CCM.O.OI  7S5329-CCM 

caul  untrminrani ,nrani ,rni 
CCS-! ,0 

IEIRN.GT.0.5!  CCS--l.fi 

TERMINAL  SPEED  rOR  INFINITE  TURm 
TSI-TS2I/SQRTI0UM1 ) 


APC  LENGTH  OF  FULL  TURN 
S-Tr-CCM 

DOES  TAHGET  ACCELERATE  or  DECEIERATE  during  turn 
FnUO«TS 1  I /TST 
IFIFNUO.GT, I ,ol  GO  TO  005 
IF ( ENUO.EG. 1 .0  I  GO  TO  006 

TIME  TO  COMPLETE  TURN.  TARGET  accelerates. 
DUMj-ATANHIFNUOI 

T7-(CR/TST).( ACOSMIC XFtS/CR I -COSH  I OUMi  |  1-DUMj  | 

Can  turn  PE  COMPLETED  RFFoRE  WEAPON  DELIVER* 

I  F ( ( TM*T 2  I «LE • T 1 1  GO  TO  00? 

course  ano  position  during  turn  at  weapon  delivery 

TaI'-(TST/CR|-(T1-TM| 

S-CR-ALOGl COSH) T  A  U ♦ 0  U  m i I/COSHIDUM) > ) 

0UM2-S/TR 

TCWD-0UH2-CCS 

XT-xT*TP-5INI0UM2» 

YT-tR-I I .P-C0SIDUM2 J )*CCS 
GC  TO  002 

SPEED  ANO  POSITION  AT  END  OF  TURN 
C07  XT-XT-TR-SINICCH) 

YT-TR-I I ,0-COSICCMI I -CCS 
TAU-ITST/CRI-T2 
TSE-TST-TANHI TAU-OUMJ I 

Reset  coefficients  tor  straight  run  aftfr  turn 

FnUO-TSF / TS2 I 

TAU-I TS21/CR) •( TI-TM-T2I 

OUM|«ATANM|rNUC» 

COURSE  AND  POSITION  AT  WEAPON  DELIVERY 
TCWO»CCH*CCS 

DUM2-TCl.  l-ALOGlCOSHlTAU-OUNl  »/CDSM(DUM|  »  ) 
XT«XT*0UM2*C0S(CCM) 

TT-TT*0UM2*S1N|CCM)-CCS 
GO  TO  0C2 


TIME  to  complete  turn,  target  speed  constant. 


227. 

006 

T?»TR.f fH/TS, I 

770. 

229. 

C 

Caw  torn  it  completed  refoRE  weapon  delivfky 

730. 

If  1  ( TM*T2  1  ,LE. T 1)  60  TO  000 

231  . 

237. 

c 

Course  »nd  position  coring  turn  at  weapon  delivery 

233. 

0Um2*CCM* l T |-TM|/T7 

23*4. 

TCWD«DUM2»CCS 

735. 

XT«XT*TP»SIN< DUh2 » 

736. 

yt»tr»m.o-cos<pum?m»ccs 

237. 

60  TO  Q02 

238. 

239. 

c 

position  at  end  or  turn 

2**0. 

008 

Xt*XT.TP«SIN(CCM) 

2H|  . 

VT»TW*l 1 .C-COSICCMI |*CCS 

292. 

293. 

c 

Reset  coefficients  for  straight  run  aftfR  turn 

29q  . 

FnU0«TSI I/TS2I 

295. 

TAD.*  TS2I/CR | « | T I-TN-T2  » 

298* 

0UMi«ATANH|FNU0» 

297. 

298. 

c 

course  ano  position  at  weapon  delivery 

299. 

TCWD«CfM»CCS 

2^0. 

0uM2»TCLI*AL06(C0SH( TAU*Du*1 I/C0SH(0UM|)> 

25|  . 

XT»xT4DUM2«C0S(CCM) 

257. 

TT*YT*DOM2*SlN(CCM).CfS 

253. 

60  TO  002 

259. 

265. 

c 

time  to  complete  turn,  target  decelerates 

756. 

005 

DuM| ■ACOTHIENUO) 

757. 

T2-(CR/TST|.(AS1NHIEXPIS/CR>»SINH|DUM1  »  ) -DUN  1  I 

258. 

2  59. 

c 

Can  turn  be  completed  BFEo»E  weapon  delivery 

260. 

IFI  (Tm*T21 .LC.TI 1  GO  TO  009 

261  . 

262. 

c 

course  and  position  during  turn  at  weapon  delivery 

263. 

TAU«tTST/CR)*tTI-TM) 

269. 

S«CR»ALOG(S|NH|  T  AU*OUm  I  )/SlNHIDllM|  1  ) 

265. 

0UM2-S/TR 

266. 

TCWD»0UM2*CCS 

267. 

XT*XT4TR»SIN( 0UM2 1 

268. 

Tt*tR*M.0-C05IDUM?||»CCS 

269. 

Go  TO  002 

270. 

271. 

c 

SPEED  and  position  AT  END  of  TURN 

27  2. 

009 

Xt-xT^TRrSINICC*!! 

273. 

VT«T**H .O-COSICCm) )«fCS 

279  . 

TaU»(TST/CR»«T2 

275. 

TsEbTST/TANHI TAU*0UN1  1 

276. 

277. 

c 

does  target  accelerate  or  decelerate  after  turn 

778. 

FnU0»TSE/TS2 I 

279. 

tr i fnuo.gt, i ,c i  go  to  oio 

280  . 

IE irNUO.EQ. 1.3)  GO  TO  01 I 

28  1  . 

2«7. 

c 

target  accelerates  after  tu**»  rfset  coefficients 

283. 

c 

Run. 
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289 
2*5 
2*6 
2*7 
2*8 
2*9 

290 
29) 

292 

293 
29*4 

295 

296 

297 
79P 
297 

300 

301 

302 
30  3 
309 

305 

306 

307 

308 


309 

310 


31  I 

312 

313 
319 

315 

316 
317 
3I« 


319 

320 

321 

322 


! 


323 

329 


325 

326 


327 

328 

329 

330 

331 

332 

333 
339 
335 


336 

337 

338 


339 

390 


TAU«<  T57I  /CR  I  •(  T  |-T*4-T2  I 
Di;M|«4TANH(fNU0) 

C  C  0  L  K  5  E  AND  POSITION  AT  WEApON  DELIVERY 

TC*D*CCM»CCS 

DuM2»TCt_l«Ai.0SlC0SH(  taU*OUM|  I/COShiDumj  i  ) 

XT*XT*PIIM2»C0S  (  C C M  ) 

Yt«TT*DUM2*S INI CCM)«CCS 
60  TO  002 

c  target  speed  constant  afte*  turn,  course  and  position  at  weapon 
C  DELJ'(ERV» 

"11  Tcwd»CC*»CCS 

DUM2»  TSC • I T|-TM-T2> 

XT«Xt*DUH2«C0S( CCM I 
Tr*rT40UM2«siN(CCM)«CCS 
GO  TO  002 

c  Tapget  decelerates  aeter  turn,  reset  coeeeicients  for  Straight 

c  Run. 

010  TaU.I TS7 I /CH ) 4 ( T 1 -TN-T2 ) 

Duf'l*ACOTM(ENUO) 

c  course  and  position  at  weapon  delivery 

TC*D«CCM4CCS 

DllM2«TCL  I  •  A  L  0  G  I  SlNH  (  TAU*DUN|  )/S|NH(DUMl  )  ) 

XT4XT4DUM2»C0S( CCH) 

YT»YT4DUM2  4S1N|CCM I  4  cc  s 

C  TeRMlNAI  for  TARGET  COURSE  AMO  POSITION 

002  continue 

C  INITIALIZE  LOOP  FOR  NUMBER  OF  SHOTS  ON  RUN 

PCHS-0.0 

Do  alt  *»1 .NSHOT 

C  WEAPON  DELIVERY  point 

Call  GA'JSS  (NR  AN  1  ,0.0,SIG0.NRAN1  ,XWD) 
call  gaussinrani,o,o,sigo,nRani,twoi 

XW0»XWD*X AP 
YwO»YwO*YAP 

c  convert  weapon  delivery  point  to  hit  probability  grid  coordinates 
OUM 1 »s I N ( TCWD ) 

Oum2«COS I TCwO I 

DuMj-XWD-XT 

OuM*t»YWO-TT 

XwDG»DUM«»4DUMi  4DUM3.DUM2 
Tw0g*DUM9  4  0UM2-DUM34DUM  j 

C  IS  THE  weapon  DELIVERY  POINT  OUTSIDE  THE  HIT  PROBABILITY  GR|D 

If  I  I XwOft.LT.GRjOXl I  I ) .OR* ( XNOG. GE • GR I D X  I  NX ) ) ,0R» 

X  (TwDG,LT,GRIDY(  |  I  ),OR,|TWDG,GE,GRlOV|NY)  I  |  GO  to  016 

c  cumulative  hit  probability  for  run 
1 1  ■  1 

020  IF  I (GR1PXI I  I ) .LE.XWDGI ,»ND, ICRIDX I  I I4|| .GT.XWOGIi  GO  TO  n 1 9 
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6 


I 


I 


1 

I 

1 


3H  . 
3*2. 
3“3. 
3MM. 
3«5. 
3M6. 
3«7. 
3H8. 
3M9. 
350. 
35\  . 

352. 

353. 
35m. 

355. 

356. 
367. 

358. 

359. 

360. 
361  . 

362. 

363. 
36M  . 

365. 

366. 

367. 

368. 

369. 


1  I » II ♦ » 

60  TO  020 
019  Jj«l 

018  IF f 1S8IDY 1 JJ| ;Lf .TWOGl ,»Nn. « GRtOT ( JJ* II .6T.TW06)  )  GO  TO  a  2 1 
J J“ J J* I 
GO  TO  0)8 

021  PCHS»PCHS*PH(  | I  , JJI-PCHS»PM(  I  I  , J 3  I 
C  LOOP  TERMINAL  FOR  HUMBER  OF  SHOTS  ON  THIS  RUN 

016  Continue 

c  cumulative  hit  prorapility  for  set 

PCH«PCH«PCHS/FNRUN 

c  loop  terminal  for  Runs  in  this  set 

017  CONTINUE 

c  print  results  for  set 

WPITEI6.22I  PCH 
WRITE16.NL2) 

WRITE! 6,76)  NX , I 6R I  OX  I t I  1 , 1  I ■ 1 ,NX I 
WPITE ( 6 ,27  1  NY, IGRIOYI I  II ,1 I»1  ,NY  > 

Nx I »NX- I 
NYIbHY-! 

WP  I  TE  (  6 ,28  1  NX1  ,NY1  ,  M  PH!  I  l  ,  JJ1  ,  I  I»|  ,HXI)  ,  JJ»|  ,NY  I) 

c  loop  terminal  for  sets  of  runs 

01 M  CONTINUE 


I 


b 


( 


1  • 

c 

2  • 

c 

3. 

c 

<*» 

c 

5» 

c 

6  • 

c 

7. 

c 

8. 

c 

9. 

c 

10. 

c 

1  1  • 

c 

12. 

c 

13. 

19. 

002 

16. 

J  4. 

17. 

COI 

18. 

19. 

?o. 

SU&RCUTINE  UNFRM| NO.NI  ,X  1  COMPUTFS  *  REAL  PSeUDO-R ANDOM  VARtABLE  X 

uniformly  distributed  over  the  interval  o.o  to  1,0  and  an  Integer 

PSEUDO-RANOON  VARIABLE  N|  UNIFORMLY  DISTRIBUTED  OVER  the  INTERVAL 
0  TO  (  ?*05-n«!  2*»8-i  1  .  IP  A  CALL  POR  t) N F R M  HAS  NOT  BEEN 
PPECEEOED  PY  ANOTHER  CALL  POP  UNFRH  OP  A  CALL  POR  GAUSS,  THp  SEED 

variable  no  must  be  an  odd  integer,  otherwise  no  should  be  the 

VALUE  OP  N l  RETURNED  "T  the  PRFCEEOiNG  CALL  FOR  UnFRM  OR  GAuSS. 
THIS  SUBROUTINE  IS  INTENDED  FOR  USE  WITH  A  COMPUTE*  HAVING  tN 
INTEGEP  WORD  LENGTH  OF  35  BITS  PLUS  SIGN  AND  AN  I NTEGER- TO-rE AL 
CONVERSION  FUNCTION  FLOATjII  WHICH  RETAINS  ONLT  THE  2 7  MS  B  f  TS 
OF  THE  INTEGER.  THE  PERIOD  OF  THE  PSEUDO  RANDOM  NUMBERS  IS 
APPROXIMATELY  2««33, 

Subroutine  unfrm<no,ni ,xi 

N I  *  I  86 36  3* No 

IFINI .6E.0)  GO  TO  001 
Nj •-) ,N| 

IF  I M ,GT .3R36973B |  I  2 )  60  TO  002 
XbFlOAT I N | I/3R3597381 T?!c 
RETURN 
ENO 


I 
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subroutine  gaussino,*m,xs,ni ,xi  computes  a  peal,  normally 
DISTRIBUTED,  PSEUDO-RANOOM  variable  X  WITH  mean  XM  and  standard 
deviation  xs,  if  a  call  foR  gauss  has  not  been  preceeded  by 
ANOTHER  call  for  GAUSS  OR  A  CALL  FOR  unprh.the  seed  variable  no 
must  BE  an  ODO  INTEGER,  OTHERWISE  no  SHOULD  be  the  VALUE  of  N| 
RETURNED  By  the  precefoing  CALL  FOR  GAUSS  OR  unfrm,  subroutine 
gauss  calls  the  subroutine  unfrm. 
subroutine  gauss<nc,*m,xs,ni  ,x » 

XbO,0 

00  OOl  I«1  ,  12 

call  unfrm<n0,no,dx > 

X  *  X  «D  X 

COI  CONTINUE 
N  i  «nO 

Xs I X-A.O I AXSAXM 

return 

end 


TERMINAL  SPEED  FOR  A  GIVEN  THRUST  LEVEL 


The  terminal  (steady-state)  speed  v  of  a  submerged  body  with  constant 
thrust  level  T  is  given  by: 


(B-1) 


where : 


p  =  density  of  the  medium 
CQ  =  drag  coefficient 
S  =  characteristic  area 


Now  let  C  denote  the  drag  coefficient  of  the  body  for  motion  along  a 
straight  path  and  c^2  the  drag  coefficient  of  the  body  in  a  turn.  Let  v^ 
denote  the  terminal  speed  of  the  body  for  a  straight  run  with  constant  thrust 
T  and  the  terminal  speed  in  an  infinite  turn  with  the  same  thrust  level  T. 
Then,  it  follows  from  equation  (B-.): 


V2  =  V1 V  C 


( B-2 ) 


Assume  CQ2  is  of  the  form: 


C  „  =  C  +  KC„ 
D2  D 1  L 


( B-3 ) 


where  K  is  some  constant  and  C  is  the  lift  coefficient  due  to  turning.  Then: 

L 


-  i  + 


(B-4) 


B-2 


The  lift  L  for  a  body  at  speed  v  during  a  turn  of  constant  radius  R  is  given 
bys 


2mv 


(B-5) 


where  m  is  the  body  mass  and  is  equal  to  the  mass  of  the  water  displaced  by 
the  body  (we  assume  trim  for  neutral  buoyancy) . 


Then,  solving  for  C  : 

Ii 


CL  -  — 
ii  pSR 


(B-6) 


Substituting  ( B-6 )  into  ( B-4 ) s 


'Dl 


1  + 


16  Km 


i-1 


'D2 


2_2_2 


( B-7  ) 


p  S  R  C. 


Dli 


Or 


Dl 


'D2 


-1 


(B-8) 


where 


16  Km 


2  2  Cn. 
p  S  Dl 


( B-9 ) 


R  is  defined  as  the  characteristic  drag  radius  squared  of  the  body.  Then, 
c 


fran  equations  (8-2)  and  (B-8): 


(B-10) 


r 

l: 


!i 


! 


1 


1 


I 


1 


Equation  (8-10)  is  used  to  compute  the  terminal  speed  for  a  target  in  a 

circular  turn  of  constant  radius  R,  at  a  constant  thrust  level  which  would 
produce  a  terminal  speed  v^  if  the  target  were  on  a  straight  run  (R  *  oo  ). 
This  formulation  is  a  consequence  of  the  fact  that  thrust  levels  for  a  ship 
or  submarine  are  ordered  in  terms  of  the  corresponding  straight-run  terminal 
speed. 

TARGET  MANEUVERS 

The  program  provides  for  a  target  maneuver  which  consists  of  an  instan¬ 
taneous  thrust  change  and/or  a  circular  turn  maneuver  with  the  start  of  the 
turn  at  the  time  of  thrust  change . 

We  define  the  following: 


v  =  instantaneous  target  speed 

v  -  initial  target  speed  at  time  maneuver  starts  ( t«0 ) 
o 

T  =  thrust  level  before  start  of  maneuver 
o 

T  *  thrust  level  after  start  of  maneuver 
T 

v  *  terminal  speed  for  the  maneuver  at  thrust  level  T 
T  T 

During  the  maneuver,  the  force  acting  on  the  target  is: 


*■  -  tt  -  !  ev2  c„2s 


«  -  '*) 

Then  the  target  acceleration  is: 


-  PC  s 
2  D2 


«  - 


(B-11) 


I 


a  *  c 


(  B- 1 2 ) 


where : 


PCD2S 

2m 


Fran  equation  ( B-8 ) : 


(B-13) 


/,  t  «C 


(B-14) 


R2 

’♦-I 

c  \  R 


The  term  1^  is  called  the  characteristic  length  of  the  target: 


2m 


pcdi! 


( B- 1 5  ) 


From  equation  (B-12)  we  obtain  the  differential  equations 


where 


cvdt  =  Kdt 
T 


and 


v 


v 


cv_ 


(B-16) 


(B-17) 


(B-18) 


The  solution  of  the  differential  equation  (B-16)  by  integration  depends  upon 
the  value  of  v.  We  identify  three  cases: 

Case  I:  Target  accelerates  during  maneuver  (v  <  1) 


Equation  (B-25)  defines  the  path  length  traveled  as  a  function  of  time.  In 
the  case  of  a  turn  maneuver,  this  is  equivalent  to  the  angle  of  turn  (course 
change),  since  the  turn  is  assumed  to  be  circular  (constant  radius). 


To  compute  the  time  required  to  travel  a  given  path  length,  or,  in  the 
case  of  a  turn,  to  complete  a  given  course  change,  we  solve  equation  (B-25) 
for  t; 


*  "i 


cosh 


■1  cs  -1 

e  cosh(tanh  \> 

c 


tanh 


v 

o 


( B-26 ) 


Target  speed  at  any  time  during  the  maneuver  is  computed  with  equation  (B-21). 
Case  II:  Target  Accelerates  during  Maneuver  (v  >  1) 

For  v  >  1,  the  differential  equation  (B-16)  has  a  solution: 


coth  ^ v  -  coth  %  =  Kt 

o 


( B-27 ) 


If  we  proceed  from  equation  (B-27)  in  the  same  manner  as  we  did  from  equation 
(B-19)  for  Case  I  (v  <  1),  we  obtain: 


1 


v  =  coth  (  Kt  +  coth  v 


( B-28 ) 


s  =  In 
c 


sinh  ( Kt  +  coth 


■’01 


-1 


sinh  (coth  v  ) 
o 


( B-29 ) 


1 

-1 

cs  /  -1  \ 

-1 

t  W  — 

K 

sinh 

e  sinh  coth  v  1 

\  °/J 

-  coth  v 

o 

(B-28) 

provides 

target  speed  at  any 

time  during 

(B-30) 


(B-29)  provides  distance  traveled  (or  course  change)  since  start  of  the  maneu¬ 
ver.  Equation  (B-30)  provides  the  time  required  to  travel  a  given  distance, 
or,  in  the  case  of  a  turn,  the  time  required  to  complete  a  given  course 


Case  III:  Target  Speed  Constant  During  Maneuver  (v  «  1) 

For  the  case  of  constant  speed  during  the  maneuver  ( \>  =  1 ) ,  the  equations 
of  motion  are  trivial,  since  we  have  linear  motion  at  constant  speed  or  motion 


along  a  circular  arc  at  constant  speed. 


APPENDIX  Cs  INPUT  DATA  VARIABLES 


NAMELIST  KL1 


NSET: 


NAMELIST  NL2 


NRUN: 


NSHOT: 


PTM: 


TS1 : 


TS2 : 


Number  of  sets  of  runs 


Number  of  runs  in  the  set 

Number  of  weapons  simultaneously  launched  for  each 
run. 

Probability  that  the  target  will  initiate  a  maneuver 
(course  and/or  thrust  change)  at  some  time  during  the 
blind  time  interval.  If  -1.0  <  PTM  <  0.0,  then  the 
target  will  initiate  a  maneuver  at  -100.0  *  PTM  per¬ 
cent  of  the  blind  time. 

Initial  target  speed  (knots). 

New  target  speed  ordered  at  time  maneuver  is  initiated 
(knots) . 


CMAX: 

CM  IN: 

TR: 

TCL: 

TR2 : 


TB: 


Maximum  magnitude  of  target  course  change  (degrees). 

Minimum  magnitude  of  target  course  change  (degrees). 

Target  turn  radius  for  course  change  (yards) 

Target  characteristic  length  (feet) 

Target  characteristic  drag  radius  squared  (yards 
squared ) 

Estimated  blind  time  (seconds) 


APR:  Aimpoint  offset  radius.  Distance  of  aimpoint  from 

estimated  target  position  at  estimated  time  of  weapon 
delivery  (yards). 

APA:  Aimpoint  offset  angle.  Relative  bearing  of  aimpoint 

from  estimated  target  position  and  course  at  estimated 
time  of  weapon  delivery,  (degrees) 

S1GTB :  Standard  deviation  of  normal  error  in  estimate  of 

blind  time,  (seconds) 

S1GTC:  Standard  deviation  of  normal  error  in  estimate  of 

target  course  ( degrees ) 

SIGTS:  Standard  deviation  of  normal  error  in  estimate  of 

target  speed  (knots) 

SIGL:  Standard  deviation  of  normal  target  localization 

errors  (yards). 

SIGD:  Standard  deviation  of  normal  weapon  delivery  errors 

( yards) . 

NX:  Number  of  hit  probability  grid  lines  perpendicular  to 

target  centerline  (program  limit:  25  maximum) 

NY:  Number  of  hit  probability  grid  lines  parallel  to 

target  centerline  (program  limit:  25  maximum) 

Seed  for  random  number  generators.  May  be  any  integer 
less  than  2^  in  magnitude. 


NRAN: 


NAMELIST  NL3 


GRIDX(I),  (1=1,  NX):  Distance  from  target,  along  target  centerline,  to  1^ 

hit  probability  grid  line  perpendicular  to  centerline. 
Positive  ahead,  negative  astern.  Array  must  be 
ordered  such  that  GRIDX(I)  <  GRIDX(I+1). 

GRIDY(I),  (1=1,  NY):  Distance  from  target  centerline  to  Ith  hit  probability 

grid  line  parallel  to  centerline.  Positive  to  port, 
negative  to  starboard.  Array  must  be  ordered  such 
that  GRIDY(I)  <  GRIDY ( 1+1 ) . 

PH ( I , J ) ,  (1=1,  NX-1),  Probability  of  hit  for  hit  probability  grid  rectangle 
(J=1,  NY-1):  defined  by  GRIDX(I)  _<  x  <  GRIDX(I+1)  and  GRIDY ( J )  <  y 

<  GRIDY  (J+1) 
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INPUT  DATA  CARDS 


D-1 


The  following  key  punch  work  sheets  represent  a  typical  set  of  input  data 
cards  for  program  execution,  with  each  line  of  each  work  sheet  representing 
one  card.  The  numbers  shown  in  columns  79-PO  are  not  part  of  the  input  data 
and  would  not  be  included  in  the  actual  input  cards.  The  numbers  are  included 
only  for  reference  purposes  in  the  text  of  this  appendix. 

Card  01  specifies  that  program  execution  will  consist  of  four  sets  of 
runs.  Parameters  for  the  initial  set  of  runs  are  contained  in  cards  02-21. 

The  hit  probability  grid  parameters,  contained  on  cards  07-21,  correspond  to 
the  hit  probability  grid  shown  on  Figure  (2). 

Input  parameters  for  the  second  set  of  runs  are  provided  by  cards  22-23 . 
These  specify  that  input  parameters  for  the  second  set  of  runs  will  be  identi¬ 
cal  to  these  of  the  first,  except  for  a  different  seed  (NRAN)  for  the  random 
number  generators.  Thus  sets  one  and  two  can  test  for  stochastic  convergence 
by  comparing  results  for  generally  different  sequences  of  random  variables 
from  the  same  statistical  distributions.  Note  that,  even  though  there  are  no 
changes  to  NL3  namelist  data  for  the  second  set,  a  blank  NL3  card  (#23)  is 
provided,  as  required,  for  proper  program  execution. 

Card  24  changes  the  magnitude  of  target  turn  maneuvers  and  the  random 
number  generator  seed  for  the  third  set.  Card  25  provides  the  required  card 
for  NL3,  even  though  there  are  no  NL3  data  changes. 

Cards  26-27  specify  that  the  input  parameters  for  the  fourth  set  will  be 
identical  to  those  of  the  third,  except  for  the  random  number  generator  seed. 
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E-3 


E-4 


ACOSH ( X ) : 


Statement  function  which  computes  the  inverse  hyper¬ 
bolic  cosine  of  X. 


ACOTH(X) : 

APAI: 

ASINH(X) : 

ATANH(X) : 

CCM: 

CCS: 

CR: 

DUM1 ,  DO M2, 

FNRUN: 

FNUO: 


Statement  function  which  computes  the  inverse  hyper¬ 
bolic  cotangent  of  X. 

Input  variable  APA  converted  to  radians. 

Statement  function  which  computes  the  inverse  hyper¬ 
bolic  sine  of  X. 

Statement  function  which  computes  the  inverse  hyper¬ 
bolic  tangent  of  X. 

Target  course  change  magnitude  ( radians ) 

Indicates  the  direction  of  the  target  course  change. 
CCM  *  +1  for  a  negative  change  (positive  y  direction 
or  left  turn).  CCM  *  -1  for  a  positive  change  (nega¬ 
tive  y  direction  or  right  turn) . 

Reciprocal  of  hydrodynamic  coefficient  c  (yards)  (see 
equation  (B-14)] 

DU  M3 ,  DUM4:  Duntny  variables  for  temporary  storage  of  intermediate 
quantities 

Floating  point  conversion  of  input  variable  NRUN 

v  -ratio  of  initial  to  terminal  speed  [see  equation 
o 

(B-20 )] 


F-2 


NRAN1 : 


Running  seed  integer  for  randan  number  subroutines 
UNFRM  and  GAUSS. 


PCHs 

PCHS: 

RN : 


S: 

SIGTCI s 

SIGTSI : 


TIs 

T2s 

I 

TAU: 

TCLI; 

I 

I 

TCWD: 

TM: 


TS1I: 

TS2Is 


l 


TSE: 


t 


I 


Cumulative  hit  probability  for  set. 

Cumulative  hit  probability  for  run. 

Randan  number  fran  a  uniform  distribution  over  the 
interval  [0,1] 

arc-length  of  a  turn  maneuver  (yards) 

Input  variable  SIGTC  converted  to  radians 

Input  variable  SIGTS  converted  to  yards  per  second 

Actual  time  to  weapon  delivery  (seconds) 

Time  required  to  complete  a  turn  maneuver  (seconds) 
Time  variable  Kt  [see  equations  (B-19)  and  (B-27)] 
Input  variable  TCI.  converted  to  yards 

Target  course  at  weapon  delivery  time  (radians) 

Time  target  maneuver  starts  (seconds) 

Input  variable  TS1  converted  to  yards  per  second 
Input  variable  TS2  converted  to  yards  per  second 

Target  speed  at  end  of  turn  (yards  per  second) 


F-3 


TST; 

XAP: 

XT: 

XTC: 

XWD: 

XWDG: 

YAP: 

YT: 

YTC: 

YWD: 

YWDG: 


Terminal  speed  for  maneuver  (yards  per  second)  [see 
equation  (B-10)] 

X- COORDINATE  OF  AIMPOINT  (YARDS) 

x-coordinate  of  target  (yards) 

Estimated  x-coordinate  of  target  (yards) 

x-coordinate  of  weapon  delivery  point  (yards) 

x-coordinate  of  weapon  delivery  point  in  hit  proba¬ 
bility  grid  coordinates  (yards) 

y-coordinate  of  aimpoint  (yards) 

y-coordinate  of  target  (yards) 

Estimated  y-coordinate  of  target  (yards) 

y-coordinate  of  weapon  delivery  point  (yards) 

y-coordinate  of  weapon  delivery  point  in  hit  proba¬ 
bility  grid  coordinates  (yards) 


